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Abstract 

It has recently been suggested that the formation of horizon-size primordial black 
hole (PBH) from pre-existing density fluctuations is effective during the cosmic 
QCD phase transition. In this Letter we discuss the dependence of PBH forma- 
tion on effective relativistic degrees of freedom, g e g during the cosmic QCD phase 
transition. Our finding is important in the light of recent cosmological arguments of 
several new classes of cosmic energy that appear from universal neutrino degeneracy, 
quintessential inflation, and dark radiation in brane world cosmology. Extra-energy 
component from the standard value in these new cosmological theories is represented 
as an effective radiation in terms of g e s . We conclude that the PBH formation dur- 
ing QCD phase transition becomes more efficient if negative extra-component of 
the cosmic energy is allowed because of the increase of the duration of the QCD 
phase transition, which leads to smaller mass scale of PBHs. This suggests larger 
probability of finding more PBHs if the dark radiation exists as allowed in the brane 
world cosmology. 
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1 Introduction 



First argument of the PBH formation was presented by Carr and Hawking (JlJ) , 
in which the key concept of the PBH formation is that the three different length 
scales of the particle horizon, the Schwarzschild radius, and the Jeans length 
of any overdense region at radiation-dominated epoch are of the same order 
of magnitude. Once an initially super horizon- size density fluctuation crosses 
the particle horizon, competition between self gravity and centrifugal pressure 
force determines whether such a fluctuation collapses into black hole or not(j2 
0). If the self gravity overcomes the pressure force, then the fluctuation will 
form a black hole with the order of horizon mass. Jedamzik (4) found that, 
during cosmic QCD phase transition, there is almost no pressure response on 
adiabatic compression of fluctuation in the mixed phase of quark-gluon plasma 
and hadron gas. This could occur because the enhanced energy density due 
to the compression of the fluctuations will turn into the high energy quark- 
gluon phase, instead of increasing centrifugal pressure force. As a result the 
overdense region from initially existing fluctuations tends to collapse efficiently 
to form the PBH during the QCD epoch when horizon mass scale is around 
1M . PBH mass spectrum is most likely dominated by the horizon mass at the 
phase transition epoch (jjj. Green and Liddle (@) verified that this speculation 
turns out to be a good approximation in their detailed studies of the PBH 
mass spectrum. 

These theoretical speculations have been studied more quantitatively in ID 
general- relativistic hydrodynamic calculations (7), where the required critical 
overdensity parameter, She = 8p/p > was found to be She ~ 0.7 for the PBH 
formation. Similar sophisticated calculations (0) have also shown numerically 
that She is even as low as 5 he ~ 0.54 by taking account of the dynamics of the 
QCD phase transition in supercooling phase by the use of EOS in the MIT 
bag model for the mixed phase. 

Recently an interesting cosmological proposal has been made that our universe 
is occupied fractionally by a new class of extra-energy component in addition 
to the ordinary matter, known as "quintessence". The quintessential scalar 
field causes the accelerating universal expansion at present (0; [IB ) and it also 
has a possibility of making significant contribution to the total energy density 



at radiation dominated epoch (11; 12; 13; 14). Such an extra-energy density 
is interpreted as to increase relativistic degrees of freedom, g eS , and thereby 
causes more rapid expansion rate at the same temperature. Similar effect is 
expected in the neutrino-degenerate universe models (flit), and the universal 
lepton asymmetry was critically studied in the context of constraining the pri- 
mordial nucleosynthesis and cosmic microwave background anisotropies. On 
the other hand, another attempt has been made quite recently to understand 
the Einstein gravity in the light of the unified theory (flU fl7l : Hit Eih . Five 
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dimensional brane world cosmology leads to an effective theory of gravity, 
which is motivated by the superstring theory or M-theory, which is a strong 
candidate for the unified theory. In this scenario the "standard" Friedmann 
equation should be slightly modified by adding several new terms which arise 
from the extra dimensions. Several authors (17; 0; Eil ) showed that there is 
a new class of effective energy among them which diminishes in proportion to 
a -4 called "dark radiation" , where a is the cosmic scale factor. The interest- 
ing point is that this dark radiation term could be even negative because it 
originates purely from fifth dimensional geometry. If this is the case, negative 
dark radiation causes the cosmic expansion slower. We can therefore interpret 
this effect as a decrease in effective degrees of freedom, g e g. 



The purpose of the present Letter is to study how strongly the change of 
the relativistic degrees of freedom g eS affects the PBH formation during QCD 
phase transition. The energy density of the universe gets larger for larger g e ^ 
at the same temperature, but the duration of QCD phase transition becomes 
shorter. These two effects operate in opposite direction to form PBHs. There- 
fore, it is not a trivial subject since many physical quantities may be affected 
in a complicated manner by this parameter g e g. In order to clarify its ef- 
fects we use a simple model for the numerical analysis instead of performing 
complicated hydrodynamic calculations. We then try to find the parameter de- 
pendence of the typical mass scale and the amount of PBHs which are formed 
during the QCD phase transition. 



2 QCD Phase Transition 



We first briefly review important quantities involved in the PBH formation 
during the QCD phase transition. Horizon mass at the QCD epoch is written 

as 

M h (T) = 0.8M© ( ) 2 (JbL-Y* (1) 

V ; ° VlOOMeV/ V51.25/ ' V ; 



where g* is the number of relativistic degrees of freedom, and T is the cosmic 
temperature. This equation shows that the horizon mass at the QCD phase 
transition is compatible to typical mass scale of MACHOs (j2oh. 



Recent progress in lattice QCD calculations and accelerator experiments have 
provided with rich information of the QCD phase transition. We employ the 
MIT bag model to express the quark-gluon plasma phase (|2ll). In this model 
the energy densities of quark-gluon plasma (p qg ) and hadron gas (ph) and the 
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pressures (p qg and ph) respectively are given by 

Pqg = ^qe T4 + 5, PhiT) = ^9hT\ (2) 
= ^p qg - 5, p h (T) = ~p h , (3) 



where g qg and are the numbers of degrees of freedom of the quark-gluon 
plasma and hadron gas phases, and B is the Bag constant. The condition of 
pressure balance between the two phases, p^Tqcd) = Viz^QCd), for a first 
order phase transition yields 



7T 2 

B = —(g qg - g h )T^ CD , (4) 
L = p qg (T QCD ) - p ft (T Q cD) = 45, (5) 



where L is the latent heat of the phase transition. This implies that the inten- 
sity of the first order QCD phase transition depends on the difference between 
the two degrees of freedom before («7 qg ) and after (g^) the phase transition. 

Duration of the phase transition is expressed by the change of the cosmic scale 
factor during the phase transition. From the conservation law of entropy we 
estimate the rate 



a 2 
a i 




1.44, (6) 



where a± and a<i are the scale factors at the beginning and the end of the QCD 
phase transition, respectively. We can learn from this that the duration of the 
phase transition is comparable to the Hubble time at the QCD epoch, and it 
is enough for fluctuations to collapse gravitationally. 



3 Analysis 



3. 1 Theoretical Model 



For the numerical analysis of the PBH formation we mainly follow the method 
developed by Cardall and Fuller (H3). The condition for the gravitational col- 
lapse of PBH is that the self gravity is stronger than the centrifugal pressure 
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force, which is expressed as 

Gp 2 S c 5 oll > P S 3 coll , (7) 



where p is energy density, p is pressure, and S'coll is the size of an overdense 
region at the turn around. This condition becomes 



Scoii ^ i / ~- ~ R.i 



Gp 



where w = - = \ in the radiation dominated era and Rj is the Jeans length. 

p 3 

Equation (8) is a well known condition for the gravitational collapse. This 
condition is also expressed, using the density contrast <5hc at horizon crossing 
(0), as 

Z w. (9) 



This means that fluctuations having density contrast above g at the horizon 
crossing can collapse to form PBH at radiation dominated era. 

We need modification of Eq.(9), taking account of the effects of the QCD 
phase transition. In order to proceed quantitative discussions, we introduce 
the quantity / and modify the above condition by 



S hc >w(l-f), (10) 

/=^A (11) 

D coll 



where S2 and Si are the sizes of initially expanding fluctuation at the beginning 
and the end of the phase transition, and S co u is the size at turn around (jiih . 
This quantity / works for softening the equation of state w and takes the value, 

< / < 1 — f-^-j • When / is the minimum value, / = 0, the fluctuation 

V ^coll J 

is not affected by the QCD phase transition at all, and Eq.(10) turns out to 
be the simplest case of Eq.(9). When / is the maximum value, f = 1 — -S 3 -, 

•^coll 

the fluctuation has been affected the most by the phase transition because 
the fluctuation starts collapsing early enough so that the overdense region 
includes the mixed phase of quark-gluon plasma and hadron gas. In the latter 
case the system does not feel any effective pressure response. S 2 and Si are 
calculated from the energy density of the background and the energy contrast 
of the fluctuation when the fluctuation crosses the particle horizon. We solve 
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the flat Friedmann equation for the background 




'da\ 2 8-kG 9 

Tt) =— (12) 



and the closed Friedmann equation for the fluctuation 



~ P S 2 -k, (13) 

- —p h K h y-g^ ~ ) ■ y M ) 

In these equations, we assume that the fluctuation has the top hat profile 
Ph = Ph(l + 8), where 5 means the density contrast at horizon crossing. We 
normalize the density evolution and fix the gauge by setting = th, Sh = Rh, 
and (dS/dr)h = (dR/dt)h when the fluctuation crosses the horizon. 

We show the quantity w(l — f) as the function of 5 in Fig.l, from which we can 
show the effect of QCD phase transition on the PBH formation. Here, x = 
Ph is the mean energy density when the fluctuation crosses the horizon, and 
pi is the mean energy density above which the universe is in pure quark- 
gluon plasma phase. This variable x indicates when the fluctuation under 
consideration crosses the particle horizon. 

In order to connect the physics of the QCD phase transition with the relativis- 
ts number of degrees of freedom, we introduce effective background degrees of 
freedom,g cff = g hg + Ag, where g^ g ~ 14.25 (7, e ± , p ± , three kinds of neutrinos 
in the standard case), and Ag is for the extra-component at the QCD phase 
transition epoch as to be discussed in the next section. The numbers of de- 
grees of freedom of the hadron phase and quark-gluon phase are respectively 
given by, = g c s + 3 (pions), and g qg = g c s + 37 (21 for quarks and 16 for 
gluons). We thus obtain the relation between the parameter x and the number 
of effective degrees of freedom g cS for horizon mass scale M, 

x _p h _ 3g c{i +U5 ( M V 2 
pi 153.75 V°- 87M ©/ ' 



3.2 Present Day PBH Density 



We consider the significance of the PBHs formed in the QCD phase transition 
in this sub-section. Present cosmological density parameter of PBHs formed 
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Fig. 1. Modification of collapsing condition for x = Qhf p\ = 15. x = 15 corresponds 
to a fluctuation with mass scale Mpbh ~ 0.2M Q (0). w{l - /), for w = \, and S 
are shown as indicated. Collapse region 6 a 5 ^ ^satisfies the collapsing condition 
of Eq.(10). 



in the QCD phase transition is represented by 
tt PBH h 2 = 4.8 x 10 6 g(T)e(T) 



T n \ 3 / T 



2.73K/ VlOOMeV 



(16) 



where h is the Hubble parameter in units of 100 km/s/Mpc, g(T) is the rel- 
ativistic degrees of freedom, and s(T) is the fraction of the radiation energy 
density which is converted into black holes at temperature T. e(T) is defined 
by 

s b 

e(T) = J F(8,T)d5, (17) 

Sa 



where F(5, T) is the probability function of a horizon volume fluctuation to 
have overdensity parameter 5 (23) 

F(5, T) = -=J- exp | ^— ] , (18) 
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/ M \ {1_n)/4 

"W = 9-5X10- t„) . (19) 



In Eq. (19) n is the spectral index of primordial density fluctuations. We 
consider the fluctuation which has overdensity larger than 5 a [see Eq. (17)] 
and can collapse into black holes, so that it satisfies Eq.(10) for a given x 
(See Fig.l). It is to be noted that even a small fraction e can result in the 
significant amount of cosmological energy density at present. From Eq.(15) we 
can estimate the mass of the PBH and the variance a for a given x and n. We 
can also calculate the efficiency e and Qpbr using Fig.l and Eqs.(15)-(19). 



4 Results and Discussions 



If the universe is modeled to have a new class of extra-energy component 
, it leads to increase relativistic degrees of freedom, g(T). In the neutrino- 
degenerate universe models (jlil ). additional degrees of freedom are written 

as 




where a's correspond to e-,/i- and r-neutrino species, and £ is the neutrino 
chemical potential divided by the neutrino temperature, £ a = p a /T a . We note 
that slightly non-zero chemical potential is favored in recent BBN analysis 
with better goodness of fit to the recent data of CMB anisotropics |l5h. which 
constrains Ag < 4.72. 



Another class of extra-energy component appears in quintessential inflationary 
model (HJ E E3; Q 



Ag = g Q = g B 



Pq 



(21) 



where ps is the background photon energy density, ps = p^o(l+ Z Y (go / g (z)) 1 ^ 3 , 
and pq is the energy density of quintessence field, pq = Qqcd/2 + ^(Qqcd)- 
Qqcd and Qqcd are the quintessence scalar field value and its time derivative 
at the QCD phase transition. Ag is constrained to be < 0.2 so as to satisfy 
the BBN and CMB analyses (Q). 



In a recent theory of brane world cosmology ( 1|| 17 ; 3 E) , completely a new 
class of cosmic energy, called dark radiation, emerges from extended Friedman 
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Fig. 2. The fraction of the radiation energy density converted into black holes in 
three cases of g e g values as a function of the mass of PBH. This function is steeply 
peaked for any g e g, which suggests that PBH forms at a particular epoch for almost 
single horizon mass corresponding to a peaked M. In this figure we have tuned 
spectral index so that ^pbh ~ 1 at present for the sake of illustration. 

equation. Its extra-component is expressed by 



Pdr 
Pb ' 



(22) 



where p D R is the dark-radiation term, p DR = ^(1 + z)\ and \i is a constant 
which comes from the electric Coulomb part of the five-dimensional Weyl 
tensor (16). We note that the dark radiation can be negative depending on 
the sign of //. This is in remarkable contrast to Eqs.(20) and (21) which are 
always positive. The BBN and CMB constraints (J2J) on Ag in the brane world 
cosmology is -0.47 < Ag < 0.12 (at the 2a C.L.). Note that only the BBN 
constraint allows —4.65 < Ag. 

We show the calculated results for g e g = 9.60, 14.25, and 20.00 in the Figs. 2 
and 3. g c s = 9.60 corresponds to p ex tra/PB = —0.27 just after the QCD phase 
transition, which is the lowest value deduced from only the BBN constraint 
in brane world cosmology (jl^). g e g = 14.25 is for the standard case, and 
g e fr = 20.00 corresponds to p e xtra/PR = 0.33 ,which is a likely case in the 
quintessential inflationary scenario (flih. Refer to ref. (flih for more critical 
discussions of the observational constraints on p ext ra/PB in quintessence scalar 
fields. 
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Fig. 3. Density parameter of PBHs OpBHat present, as a function of the spectral 
index n of the primordial density fluctuations. 

The key results are summarized as follows. Fine tuned Gaussian blue spectra 
fluctuation is needed with spectral index n « 1.354 ± 0.005, 1.377 ± 0.005, 
1.395 ±0.005 in the cases of g e s = 9.60 , 14.25, and 20.00 respectively, in order 
to have significant PBHs energy density Opbh at present, i.e. 0.01 < f2pBH < 1- 



Typical mass scales of PBH are O.116M , O.136M , and O.159M for g e g = 
9.60, 14.25 and 20.00, respectively, where three mass scales refer to maximum 
e(T)'s. (Eq.(16) and Fig.2). The spectral indices as shown in Fig. 3 are com- 
patible with those inferred from the data of COBE observations, but larger 
than the best fit analysis of the combined different cosmological data sets 
(J2^). These indices are slightly larger than the values constrained from the 
evaporation of PBHs by Hawking radiation in much smaller mass scale (|23h. 
Increasing the effective number of degrees of freedom g e g causes the shorter 
duration of the QCD phase transition because the universal expansion obeys 



o-i \gh) V 9cs + 3 J 



(23) 



instead of Eq.(6). Accordingly, the duration of the mixed phase also becomes 
shorter, making less efficient PBH formation (Fig.3). The typical mass scale of 
the PBHs, however, increases with increasing g e d (Fig.2). This is because the 
fluctuations that are most effectively influenced by the QCD phase transition 
enter the horizon at later time with increasing number of degrees of freedom. 
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5 Conclusion 



We studied the formation mechanism of PBHs during the QCD phase tran- 
sition, being motivated by the cosmological interest in MACHOs as a dark 
matter and also new classes of cosmic extra-energy components for the de- 
generate neutrinos, dark energy, and dark radiation. In the present paper we 
take account of the effects of these extra-components by means of parame- 
terizing their effect as an increase or decrease in effective relativistic degrees 
of freedom. We conclude that the PBH formation during QCD phase tran- 
sition becomes more efficient with an negative extra-component because of 
the increase of the duration of the QCD phase transition. We also found that 
in such a case the typical mass scale of PBH becomes smaller. These results 
suggest larger probability of finding more PBHs that have typically sub-solar 
mass M ~ 0.1 M Q if the dark radiation exists as allowed in the brane world 
cosmology. 
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